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Abstract 

A symmetry reduction of the Dirac equation is shown to yield 
the system of ordinary differential equations whose integrability by 
quadratures is closely connected to the stationary mKdV hierarchy. 

Consider the Dirac equation of an electron 

»E 7 M VV- (e£ + il> = 0, (1) 

moving in the electric field 

A = A (x 3 ), A 1= A 2 = A 3 = 0. 

In the above formulae 7 M are 4x4 Dirac matrices, ip = i/j(xq, Xi,x%, x 3 ) is a 
four- component complex-valued function and e, m are constants. 

The form of the vector-potential A^ imply the following Ansatz for the 
spinor field ijj(x): 

ijj(x) = (f(x 3 ). 

Inserting this expression into the Dirac equation ([!]) yields system of ordinary 
differential equations (ODEs) for the four-component function f(x 3 ) 



(ff - (2e7 3 7oA) - imj 3 )(p = 0. 
* e-mail : rzhdanov@apmat . freenet . kiev.ua 



(2) 



Denoting 



x = 2x 3 , V(x) = eA (x 3 ), 

1 i i 

Ji = 2^0, J 2 = ^73, J? = 2^370, 

we rewrite (0) in the following form: 

Cip = (D x -V(x)J 3 -mJ 2 )ip = 0, D x = ^-. (3) 

Note that the 4x4 matrices J\, J2, J3 fulfill the commutation relations of 
the Lie algebra so(3) 

[J a , J b ] = J c , (a, b, c) = cycle (1, 2, 3). (4) 

To integrate system of ODEs (0) we will make use of its symmetry prop- 
erties. The general routine for calculating symmetry group admitted by a 
differential equation is the infinitesimal Lie method. This method makes it 
possible to reduce the problem of constructing the maximal symmetry group 
to integrating some linear system of partial differential equations (called de- 
termining equations). The general solution of the latter gives rise to the 
maximal transformation group admitted by the equation under study (for 
more detail, see, e.g. (0,0). However direct application of the Lie method to 
system (|3|) is inefficient, since integration of the corresponding determining 
equations is, in fact, equivalent to integration of the initial system of ODEs 
©• 

That is why to be able to integrate the determining equations one has 
inevitably to impose some a priori restrictions on the choice of symmetry 
operators. Our idea is to look for a symmetry operator Q admitted by system 
(|3p in the form of an nth order polynomial in m with matrix coefficients. 
First, we consider in some detail the case when n = 3 and then give the 
results on calculating Lie symmetries obtained for the case of an arbitrary 
nGN. 

Thus, we adopt for a symmetry operator the following Ansatz: 
3 

Q = Y1 ( a k( x ) + h(x) m + c k (x) m 2 + d k (x) m 3 )J k , (5) 
fc=i 

where a k , b k , c k , d k are some smooth complex- valued functions. 
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Inserting the expression for Q into the invariance criterion [C, Q] — and 
splitting with respect to the powers of m and then with respect to linearly 
independent matrices J\, J 2 , J3 we get the system of determining equations 
for the functions a/., b}., Ck, d& 

d\ = 0, d$ = 0, 

4 - v% = 0, 4 + ra 2 - c 3 = 0, 4 - ci = 0, 

c' 2 - 7ci = 0, c[ + Vc 2 - 63 = 0, 4 - 61 = 0, 
- = 0, + y& 2 - a 3 = 0, 4 - 01 = 0, 

4 — Vai = 0, 4 + ^ a 2 = 0, 0,3 = 0. 

Integrating the above system of ODEs yields 

d± = 0, d 2 = C\, ds = 0, 
ci = 0, c 2 = C 2 , c 3 = (7iV, 

61 = W, 6 2 = l -c x v 2 + C 3 , 63 = C 2 V, 

ai = C 2 V, a 2 = X -C 2 V 2 + C 4 , a 3 = + iy 3 ) + C 3 V, 

where C\, C 2 , C3, C4 are arbitrary constants and furthermore the potential 
V(x) has to satisfy the following nonlinear ODEs 



Thus we have established that if the function V(x) is a solution of the 
stationary mKdV equation 



This symmetry solves the problem of integrability of system of ODEs (H) 
by quadratures due to the assertion given below. 





(6) 



then the initial system of ODEs (§) admits the Lie symmetry 
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Lemma 1 Let the system of ODEs 

Ci;=(J- + f a {x)Q^ j = 0, (7) 

where Qi,Q2,Q3 are constant matrices forming a basis of the Lie algebra 
so(3), admit a Lie symmetry 

3 

x = J2 9a(x)Q a . 

a=l 

Then it is integrable by quadratures. 

Proof. Making a change of dependent variables 

ifj -> $ = V(x)lp, V(x) = exp j^j ^ a (^)Qa| 

we can always transform the operator X to become 

X = V- 1 XV = g(x)Q 1 , g(x)^0 

and what is more this transformation preserves the structure of system ([?]). 
The invariance criterion [C, X] = 0, where 

I = v~ l CV = £ f a (x)Q a , 

0=1 

implies that 

g(hQz - / 3 Q 2 ) - g'Qi = 0. 

As the matrices Qi,Q2, Q3 are linearly independent, hence it follows that 

h = °) /s = 0, £ = const. 
Consequently, the transformed system of ODEs necessarily takes the form 

and is evidently integrable by quadratures. The lemma is proved. 
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Hence, we get a remarkable fact: if V(x) is a solution of the stationary 
mKdV equation (j^), then system of ODEs (fjj is integrable by quadratures. 

Now we turn to the case when a Lie symmetry is looked for as a polyno- 
mial in m of an arbitrary order n 

Q=i:i:f k a {x)j a m n - k . 

k=0 a=l 

The invariance criterion [C, Q] — yields the following system of deter- 
mining equations for the coefficients of the operator Q: 

/? = 0, /? = 0, 

(/ 3 fc )' - fi +1 = o, U'i)' - vf k = o, 

(/* )' + Vfi - f! +1 = 0, fc = 0,...,n-l J 
(/o)' = 0, {K)'-Vft = 0, {ffi' + Vft = Q. 

We have obtained two classes of solutions of the above system of ODEs 
which are given below 

I. n = 2N + l, Ne~N, 
/o = 0, / 2 = l,/ 3 ° = 0, 

f? +1 = /| fc+1 = 0, f* +l = R k , fc = l,...,JV, 
f 2k+2 = D x R k , fl k+1 = {V- £>~ V) R k , 



13 

and the equation 



/| fc+2 = 0, fc = 0,...,JV-l 



D x Rn = (8) 



holds. 



n = 27V + 2, iVGN 
/o = 0, / 2 = l,/ 3 ° = 0, 

= f? +1 = 0, /| +1 = ^, fc = l,...,JV 
f 2k+2 = D x R k , fl k+x = (V - D~ V) ifc, 



'3 

and the equation 



i2 



JV+l 



5 



holds. 

In the above formulae we make use of the following notations 

k 

R k = J2 C 3 (D 2 X + V 2 -VD~ l Vy V, k = 0, . . .N + 1, 

3=0 

where C , . . . , C^+i are arbitrary real constants and D~ l is the inverse of 
D x . 

A reader familiar with the soliton theory will immediately recognize the 
operator X = D 2 + V 2 — V D~ 1 V X as the recursion operator for the mKdV 
equation ||, [|] 

v t + V xxx + ^V 2 V X = 0. 

Acting repeatedly with the operator X on the trivial conserved density Iq = 
V we get the whole set of conserved densities of the mKdV equation. Next, 
the operator 

y = D X XD~ X = D\ + V 2 + V X D~ 1 V 

is the second recursion operator for the mKdV equation. Its repeated action 
on the trivial Lie symmetry Sq = V x yields the whole hierarchy of the higher 
symmetries of the mKdV equation. Hence it follows, in particular, that the 
condition (|8|) is rewritten in the form 

N 

Y,c k s k = o, s k = y k v. (9) 

k=0 

The above equation is nothing else than the higher stationary mKdV 
equation. Provided iV = 1 it reduces to the standard stationary mKdV 
equation (^|). 

Hence, due to Lemma 1 it follows the validity of the following assertion. 

Theorem 1 Let the function V(x) satisfy the higher stationary mKdV equa- 
tion (||) with some fixed N and Cq, . . . , C^. Then, the system of ODEs (^) 
is integrable by quadratures. 

It is a common knowledge that the stationary mKdV hierarchy is reduced 
to the stationary KdV hierarchy with the help of the Miura transformation 
(see, e.g. [[|, ||). Furthermore, the latter are integrated in terms of ^-functions 
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[|]. Consequently, the system of ODEs (0) is also integrable by quadratures 
thus giving rise to exact solutions of the initial Dirac equation ([!]). 

Thus symmetry analysis of a very simple reduction of the Dirac equation 
(there are quite a few of much more sophisticated reductions, see 0) reveals 
such important elements of the inverse scattering technique for the stationary 
mKdV equation as the recursion operators, infinite number of conserved 
densities, the hierarchy of higher symmetries of the mKdV equation which 
makes the integration of system (0) very rich in results. For those involved 
into the inverse scattering business this is not surprising at all, since the above 
used procedure is just an inversion of the famous approach to analysis of 
solitonic equations via the Lax pair || . More exactly, we reduce the problem 
of integrating linear equation to solving the nonlinear one. The crucial point 
is that the this nonlinear differential equation can be integrated which enables 
us to construct the general solution of the initial linear equation. 

Now it seems reasonable to carry out systematic analysis of all reductions 
of the Dirac equation by inequivalent subgroups of its symmetry group. We 
believe that in this way other hierarchies of higher solitonic equations will be 
obtained. This situation is in some analogy to symmetry reductions of the 
self-dual Yang-Mills equations yielding almost all known integrable solitonic 
equations (see, e.g. §|). 

One more important point is that the Lie transformation groups gen- 
erated by symmetries Q obtained above are not subgroups of the maximal 
symmetry group C(l, 3) (g> U(l) admitted by initial system ([[]). These sym- 
metries correspond to conditional symmetry of the Dirac equation (Jl]). Let us 
note that conditional symmetry of the linear and nonlinear Dirac equations 
was studied in J7L llOL pj|. 
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